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Abstract: Optical coherence elastography employs optical coherence 
tomography (OCT) to measure the displacement of tissues under load and, 
thus, maps the resulting strain into an image, known as an elastogram. We 
present a new improved method to measure vibration amplitude in dynamic 
optical coherence elastography. The tissue vibration amplitude caused by 
sinusoidal loading is measured from the spread of the Doppler spectrum, 
which is extracted using joint spectral and time domain signal processing. 
At low OCT signal-to-noise ratio (SNR), the method provides more 
accurate vibration amplitude measurements than the currently used phase- 
sensitive method. For measurements performed on a mirror at OCT SNR = 
5 dB, our method introduces <3% error, compared to >20% using the phase- 
sensitive method. We present elastograms of a tissue-mimicking phantom 
and excised porcine tissue that demonstrate improvements, including a 50% 
increase in the depth range of reliable vibration amplitude measurement. 
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1. Introduction 

Optical coherence elastography (OCE) is a technique for imaging the microscopic mechanical 
properties of tissue [1], with the potential to provide at least an order of magnitude 
improvement in spatial resolution over elastography based on ultrasound [2] or MRI [3]. In 
OCE, a sample is subjected to mechanical loading (stress) and the resulting local displacement 
is measured using optical coherence tomography (OCT) [4]. The strain, given by the slope of 
the displacement versus depth, is displayed in a two- or three-dimensional elastogram. Errors 
in the displacement measurement are amplified in the strain calculation [5], degrading the 
elastogram. New methods to measure displacement with improved accuracy and dynamic 
range will improve elastogram quality and increase OCE's potential for clinical utility. 

OCE may employ either quasi-static (<10 Hz) or dynamic (>10 Hz) loading. Initially, 
quasi-static loading was employed and the displacement between consecutive B-scans was 
measured using cross correlation-based speckle -tracking algorithms [1,6,7]. However, the 
dynamic range of such methods is limited by speckle decorrelation [8-10]. Quasi-static, 
phase-sensitive OCE was proposed, in part, to overcome this limitation by improving the 
lower limit on displacement [11]. This method is based on the linear relationship between the 
phase of the OCT signal and sample motion [12,13]. Recently, a needle -based extension to 
quasi-static, phase-sensitive OCE has been reported [14]. 

In dynamic OCE, quasi-static loading and displacement measurement are replaced by 
sinusoidal loading and time-dependent displacement (vibration amplitude) measurement. 
Dynamic OCE has also employed a phase-sensitive method [15-18], which is similar to 
methods used in magneto-motive [19] and photothermal OCT [20], and to measure vibration 
in the ear [21]. For in vivo imaging, dynamic methods provide the opportunity to reduce 
motion artifacts by vibrating the tissue at frequencies above those caused by subject motion 
[18,22]. Dynamic phase-sensitive OCE also offers several alternative modes of measurement. 
Variants include: spectroscopic OCE that generates contrast through probing changes in tissue 
mechanical resonances [16]; measurement of shear wave velocity to extract elasticity [23-25]; 
and internal loading using acoustic radiation force [26,27]. 

Despite the prevalence of the phase-sensitive method in the dynamic OCE techniques 
reported to date, it has two limitations. Firstly, because the minimum detectable phase 
difference, i.e., the phase sensitivity, is inversely proportional to the signal-to-noise ratio 
(SNR) of the OCT signal [12], the displacement precision reduces with depth in a scattering 
sample. This dependence limits the range of precise measured strain to significantly less than 
the OCT depth range. Secondly, the decreasing phase sensitivity with depth causes the 
incidence of phase wrapping to markedly increase [28]. Phase wrapping causes the mean 
measured displacement to tend towards zero, which leads to an underestimation of the true 
value [29]. This systematic error severely limits elastogram accuracy and quality. 

In this paper, we propose a new method for improved vibration amplitude measurement in 
dynamic OCE, which is based on an analysis of the frequency tones of the Doppler spectrum 
obtained from a sample subjected to sinusoidal loading. Since it is related to joint spectral and 
time domain (STd) OCT, previously proposed for flow imaging [29,30], we refer to it as joint 
spectral and time domain optical coherence elastography (STdOCE). STdOCE provides more 
accurate vibration amplitude measurement in the case of low OCT SNR (<20 dB), which 
extends the depth range and significantly improves the elastogram quality. Elastograms of 
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excised porcine tissue demonstrate a 50% increase in the depth to which strain is accurately 
measured. 

The paper is organized as follows. In Section 2, we describe the theory of vibration 
amplitude and strain measurement for STdOCE and phase-sensitive OCE. In Section 3, we 
present the experimental setup and method. In Section 4, we demonstrate the advantage of 
STdOCE using results from a vibrating mirror, a tissue-mimicking phantom containing a stiff 
inclusion, and a section of freshly excised porcine tissue. In Sections 5 and 6, respectively, we 
present a discussion of the results and the main conclusions. 



In STdOCE, two Fourier transforms are performed on data acquired with a spectral-domain 
(SD) OCT system. The first Fourier transform is performed in the wavenumber domain to 
retrieve the depth-resolved OCT signal. The second is performed in the time domain and, for a 
sample voxel containing vibrating components, results in a Doppler spectrum containing 
frequency tones. In STdOCE, these frequency tones, described by Bessel functions of the first 
kind, are used to extract the vibration amplitude at each depth in the sample. 

STdOCE is related to a previously reported dynamic OCE method based on time -domain 
(TD) OCT [22,31]. That method uses the ratio of odd and/or even frequency tones to extract 
vibration amplitude. The bandwidth of these tones in TD-OCT increases with both the 
bandwidth of the optical source and the A-scan velocity [31,32], which causes them to overlap 
in frequency and restricts measurements to only the first few tones. This restriction limits the 
accuracy of the method, and requires the use of a prohibitively slow A-scan acquisition rate 
(~1 Hz). Furthermore, difficulties arise when a ratio of frequency tones approaches an 
asymptote caused by a vanishing frequency tone, unless complex post-processing techniques 
are employed [33]. STdOCE overcomes these limitations. The increased number of tones in 
the Doppler spectrum provided by SD-OCT enables use of a method closely related to one 
proposed in ultrasound elastography [34], in which the vibration amplitude is determined from 
the variance or spectral spread of frequency tones in the Doppler spectrum. 

2. 1. OCT signal of a vibrating sample 

Consider a spectral-domain OCT system consisting of a stationary mirror in the reference arm 
and a mirror in the sample arm oscillating sinusoidally in the axial direction. After performing 
a Fourier transform on the detected OCT signal in £-space [4], its amplitude may be expressed 
as [31] 



where A z is proportional to the amplitude in the sample and reference arms and the argument 
of the cosine function, <p toi , represents the total phase of the signal. cp DC is the quasi-static 
phase of the interferometer and f s is the vibration frequency of the mirror. The vibration 
amplitude, d s , can be determined from the maximum phase deviation, <p s , according to [3 1] 



where 6 is the angle between the incident beam and the vibrating mirror and A 0 is the mean 
wavelength of the light source. The displacement of the mirror may be described by 



2. Theory 



s(z, t) = A z cos [<p M (z, 0] = A z cos [<p DC + <p s (z, t) sin (2xfj)] , 



(1) 



4n;d s (z,t)cos0 



(2) 



d(z,t) = d s (z,t)sm(2xf s t), 



(3) 



and the corresponding velocity may be expressed as 
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v(z,0 = 2*-/,</,(z,0cos(2*-/,f). (4) 
The Doppler shift introduced to the detected OCT signal is given by [4] 

2v(z^)cos6> 

Substituting Eq. (4) into Eq. (5) allows the Doppler shift to be written in terms of the vibration 
amplitude, d s , as 



f D ( Z ,t) - f m OOB(^) = 4 »/.W)CO.g cQ<2nfstl 



(6) 



where f m is the maximum Doppler shift, which is related to the phase defined in Eq. (1) by 
[35,36] 

f ( t \ 1 d( Ptot (n. 

Equations (6) and (7) can be used to rewrite the argument of the cosine function in Eq. (1) in 
terms of the maximum Doppler shift: 

s(z,t) = A z cos[^ DC + /?sin(27r/ s ,?)], (8) 

where /?is the modulation index, defined as f m I f s [35,36]. From Eq. (8), we observe that the 
amplitude of the OCT signal is modulated by the Doppler shift introduced by the vibrating 
mirror. The vibration amplitude of the mirror can be determined from a measurement of p. To 
demodulate the signal, i.e., extract p from s(z,t) , we employ an algorithm closely related to 
one previously presented in the related field of sonoelasticity imaging [34]. 

2.2. Measuring vibration amplitude from the Doppler spectrum 

Using trigonometric identities, it is possible to express Eq. (8) as a series of Bessel functions 
of the first kind [34,35]: 

co 

s( Z ,t) = A z y £j n (j3)cos(n2xf s t), (9) 

where J„ represents the «th-order Bessel function. Importantly, Eq. (9) shows, in the case of a 
sinusoidally vibrating mirror, that the amplitude is not only modulated at the fundamental 
vibration frequency, f s , but also at the overtones, nf s , each weighted by the corresponding 
Bessel function, J n . Fourier transforming Eq. (9) allows us to determine the Doppler power 
spectrum, given by 

A 2 00 

S(z,f) = ^- y £j 2 n (/3)[S(f-nf s ) + 6(f + nf s )]. (10) 

As described previously [34], the vibration amplitude is estimated from the positive half- 
spectrum. The theoretical temporal waveform and positive spectral components of the 
vibrating signal in a single voxel are plotted in Fig. 1 based on a vibration frequency of 50 Hz 
and three values of p: 1.5, 3 and 6 (shown increasing from top to bottom in the figure). 

For the spectrum described by Eq. (10), the mean frequency is zero, since 
J-niP) ~ J n (P) > an d the power spectrum is, thus, symmetric about zero frequency. The 
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Fig. 1 . Theoretical ((a)-(c)) temporal and ((d)-(f)) spectral components of an OCT signal from a 
miiTor vibrating at 50 Hz with P equal to (a), (d) 1.5; (b), (e) 3; and (c), (f) 6. 

spectral spread (or variance), aj- , may be calculated from the zeroth and second moments of 
the spectrum, m 0 and m 2 , respectively, according to [34,37] 



(11) 



For sinusoidal modulation, it has been shown previously [34] that 

n— co 

«o = 2>„ 2 (/?)=l, 

n-— ao 

and 



in 



= Z(«/j 2 >/ n 2 (A) = 



ififf 



(12) 



(13) 



Substituting Eqs. (12) and (13) into Eq. (11) allows us to express the modulation index, /?, in 
terms of the spectral spread [34]: 



a/2ct 



/ 



(14) 



Having determined /? from the Doppler spectrum, the vibration amplitude of the sample 
can be calculated using the analysis presented in Section 2.1. Until now, to simplify the 
analysis, a vibrating mirror has been considered. However, this analysis is also applicable to 
turbid samples under sinusoidal loading. In this case, d t is measured at each depth in the 

sample and the strain, e, may be calculated from [22] 
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where Ac/ s is the change in vibration amplitude over an axial range Az within the sample. In 

this paper, we calculate the strain using a weighted least squares strain estimation method 
described previously [5]. We then generate 2D elastograms by scanning the beam across the 
sample in the x-direction. 

2.3. Phase-sensitive vibration amplitude measurement 

To allow us to compare and contrast STdOCE with phase-sensitive OCE, in this subsection 
we briefly describe the phase-sensitive method of measuring vibration amplitude, previously 
presented in [21]. Consider again the case of a sinusoidally oscillating mirror described in 
Subsection 2.1. From Eqs. (1) and (2), the total phase of the detected OCT signal is given by 

4xd f (z,t)cos8 / \ 
<P M = <Pdc+ ^r 1 srn(2^//). (16) 

Differentiation can be performed to eliminate the phase term cp DC , giving 

. , , %x 2 d(z,t)fAtcos0 ,„ . . 
A^(z,0 = cos {2af,t), (17) 

where At is the A-scan acquisition time. From Eq. (17), the phase difference between 
successive A-scans can be used to extract the vibration amplitude, d t (z, t) . It has been 

demonstrated that the sensitivity of this method can be improved by performing a Fourier 
transform on the phase difference, which gives [21] 

p(z, f) = FT[(A^(z, 0)] = d s (z, /, )/, At cos 0 [S(f - f s ) + S(f + /,)] , (1 8) 

where p(z,f) is the power spectral density of the phase signal. The sample vibration 
amplitude at the known frequency, f s , can be calculated from the following expression: 

(2x)f s Atcos0 

If d s is measured at each depth in the sample, Eq. (15) can be used to calculate the strain. 
It is important to note that Eq. (19) is valid if the phase difference is correctly unwrapped. For 
significant wrapping errors, caused by the presence of noise, d s is underestimated. In the next 

section, we describe the experimental setup and methodology used to implement the analyses 
described in this section. 

3. Methods 

3. 1. Experimental setup and procedure 

Results were obtained using a spectral-domain OCE system employing a sup erlumine scent 
diode source with mean wavelength of 835 nm and 55.4 nm bandwidth, illuminating the 
sample with 10 mW of optical power. The theoretical lateral resolution is 15 urn and the 
sensitivity was measured to be 105 dB for an exposure time of 50 (is. The spectrometer 
consisted of a diffraction grating with 1800 lines/mm, a lens with a focal length of 50 mm and 
a line-scan camera with 4096 pixels. The sensitivity roll-off was measured to be 4 dB/mm. 
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The camera exposure time was 50 |is and the A-scan acquisition time was 500 |is. Lateral 
scanning over a 3 mm range was performed using a galvanometer mirror, driven by a 
sawtooth waveform. The phase sensitivity of the system is 9 mrad, for OCT SNR = 50 dB, 
calculated from the standard deviation of phase measured from a stationary mirror [12]. This 
value exceeds the theoretical minimum of 2.2 mrad because of galvanometer jitter. 

Samples were mechanically loaded using a piezoelectric ring actuator, described 
previously [22], that enabled imaging and compression from the same side of the sample 
(from the top in all presented elastograms). For each sample, a refractive index of 1.4 was 
assumed. Samples were preloaded by compression against a rigid upper boundary to ensure 
full uniform contact. This was achieved by displacing the compression plate by 250 urn 
beyond the point of first contact. Each sample was tilted with respect to the beam axis to 
minimize specular reflection from the air-sample interface. To maintain uniform compression, 
the ring actuator compression plate and the upper rigid boundary were tilted by the same angle 
in both x and y planes. In all cases, a 50 Hz sinusoidal wave was applied to the actuator, 
resulting in unloaded vibration amplitudes of up to 0.45 urn. 

3.2. Phantom fabrication 

Tissue-mimicking phantoms with dimensions (length x width x height) of 10 mm x 10 mm x 
2.75 mm were fabricated using room-temperature vulcanizing silicone as the bulk medium 
and titanium dioxide (Ti0 2 ) as optical scatterers [38]. Two phantom designs were fabricated. 
Phantom 1 was designed to be optically and mechanically homogeneous. It was fabricated 
from soft silicone (Elastosil® P7676, Wacker, Germany) and contained 1.5 mg/ml of Ti0 2 
particles. Its elastic modulus was measured using an Instron materials testing system to be 25 
kPa. Phantom 2 was designed to be a soft bulk medium containing hard inclusions. It was 
fabricated from the same silicone used in Phantom 1. Inclusions of measured average 
diameter of 0.4 mm were fabricated using the method presented in [5]. The inclusions 
contained 2.5 mg/ml Ti0 2 particles and were fabricated using a harder silicone (Elastosil® RT 
601, Wacker, Germany), with measured elastic modulus of 4 MPa. They were embedded in 
the soft surrounding silicone as described previously [5]. The base of an inclusion was located 
~0.5 mm below the phantom surface. 

3.3. Data Processing 
3.3.1. STdOCE 

Figure 2 shows the steps in the data acquisition and processing procedure. Phantom 1 was 
placed under load as described in Subsection 3.1. In Fig. 2(a), 2000 spectral interferograms 
acquired from the same lateral position in Phantom 1 are shown. These spectra are modulated 
in both wavenumber and time: modulation in wavenumber is caused by scattering from the 
sample at each depth, whilst modulation in time is caused by loading-induced motion, visible 
in the inset of Fig. 2(a). In Fig. 2(b), A-scans obtained by performing a Fourier transform in 
wavenumber are shown. Each A-scan corresponds to the standard Fourier transform 
performed in SD-OCT to obtain depth-resolved measurements [4]. The sample motion at the 
position indicated by the green box is visible in the inset of Fig. 2(b). 

In STdOCE, a second Fourier transform is performed in the time domain, i.e., at each 
depth in the set of A-scans, in order to access the Doppler spectrum. In Fig. 3(c), the Doppler 
spectrum due to motion at the depth indicated by the red line in Fig. 3(b) is shown, with nine 
overtones in evidence. The frequency resolution of the sampled Doppler spectrum is 
determined by the ratio between the A-scan sampling frequency and the number of points, i.e., 
A-scans, in the Fourier transform. Furthermore, the highest measurable frequency is set by the 
Nyquist sampling criterion, i.e., half the A-scan sampling frequency [39]. Also, the range over 
which the frequency tones can be unambiguously detected is halved by the well-known 
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Fig. 2. Data processing procedure for STdOCE: (a) 2000 spectral interferograms resampled to 
ir-space; (b) A-scans obtained after Fourier transform in £-space; (c) Doppler spectrum after 
Fourier transform in time, performed at depth indicated by the red line in (b); and (d) Vibration 
amplitude calculated using the spectral spread algorithm. The red dot in (d) corresponds to the 
vibration amplitude calculated from the Doppler spectrum in (c). 

complex ambiguity problem [40]. These points need to be considered to ensure the highest 
overtone is adequately sampled. 

Having calculated the Doppler spectrum, the spectral spread algorithm described in 
Subsection 2.2 is used to calculate the vibration amplitude. Specifically, a f is calculated 

from the Doppler spectrum and then substituted into Eq. (14) to calculate /?. The vibration 
amplitude calculated from the spectrum in Fig. 2(c) is indicated by the red dot in Fig. 2(d). 
This procedure is repeated for each depth in the sample, allowing a vibration amplitude plot 
(blue line in Fig. 2(d)) to be generated. As expected for this mechanically homogeneous 
sample, the vibration amplitude decreases approximately linearly with depth. 2D vibration 
amplitude images and elastograms can be obtained by repeating this procedure whilst 
scanning the beam over the surface of the sample, which is demonstrated in Section 4. 

3.3.2. Phase-sensitive OCE 

In Fig. 3, the same data used to generate the results shown in Fig. 2 are used to describe the 
processing procedure for phase-sensitive OCE. Figure 2(a) is reproduced in Fig. 3(a), for 
convenience. Figure 3(b) shows the image generated by calculating the phase difference 
between consecutive A-scans, after performing a Fourier transform in wavenumber. The 
bright/dark banding in time is caused by loading-induced motion. Figure 3(c) shows the 
spectrum of the phase difference obtained by performing a Fourier transform in time, at the 
depth indicated by the red line in Fig. 3(b). The red lines in Fig. 3(b) and Fig. 2(b) indicate the 
same depth. In this spectrum, the peak corresponding to the vibration frequency is clearly 
visible at 50 Hz. A second peak, due to galvanometer jitter, is also visible at approximately 
160 Hz. The vibration amplitude is calculated from the spectrum using Eq. (19). The vibration 
amplitude corresponding to the spectrum in Fig. 3(c) is indicated by the red dot in Fig. 3(d). 

This procedure is repeated for each depth in the sample, allowing a vibration amplitude 
plot (blue line in Fig. 3(d)) to be generated. Similarly to Fig. 2(d), in Fig. 3(d) the vibration 
amplitude decreases approximately linearly with depth from an initial value of -0.35 urn. 
Comparing the two figures, it can be seen that the accuracy of the vibration amplitude 
measured using the phase difference method is much lower than the STd method at depths 
>0.8 mm, as indicated by large negative fluctuations in vibration amplitude in Fig. 3(d). This 
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Fig. 3. Data processing procedure for phase-sensitive OCE: (a) 2000 spectral interferograms 
resampled to A-space; (b) Phase difference image obtained after Fourier transform in A-space; 
(c) Phase difference spectrum after Fourier transform in time, at depth indicated by the red line 
in (b); and (d) Vibration amplitude calculated using the analysis presented in [21]. The red dot 
corresponds to the vibration amplitude calculated from the phase difference spectrum in (c). 

is related to the probability density functions (PDFs) for phase and amplitude and is discussed 
in detail in Section 4.1. 

3.3.3. Elastogram generation 

Elastograms were generated from 128 vibration amplitude plots, recorded whilst scanning the 
beam over a lateral range of 3 mm. Each vibration amplitude plot was generated using 1 000 
A-scans. To generate elastograms, the strain was calculated over an axial range of 75 urn. 
Each pixel in an elastogram corresponds to the strain calculated at that location normalized to 
the maximum strain in the elastogram. 

4. Results 

4.1. Vibrating mirror 

Figure 4 shows the measured vibration amplitude of an oscillating mirror, determined using 
the STd and phase -sensitive methods on the same 2000 A-scans, versus the OCT SNR, which 
was varied using neutral density filters. Figure 4 shows the mean vibration amplitude (±1 
standard deviation) calculated from ten consecutive measurements. For the STd method, we 
observe a relatively constant vibration amplitude; whereas, the phase -sensitive method 
significantly underestimates the vibration amplitude for low OCT SNR. At the highest OCT 
SNR of 55 dB, the mean vibration amplitude of the mirror measured using both methods was 
0.43 um. The standard deviation of the phase-sensitive and STd methods was 1.5 nm and 2.7 
nm, respectively. As the standard deviation may be considered to represent the vibration 
amplitude sensitivity [12], this result indicates that for high SNR, phase-sensitive OCE 
provides more precise measurements. This result is explained by considering the respective 
PDFs of phase and amplitude [29]. For high OCT SNR, the phase PDF tends towards a delta 
function, whilst the amplitude PDF is approximated by a Gaussian distribution, with variance 
determined by the shot noise of the system. However, the situation changes dramatically as 
the OCT SNR decreases. The phase PDF tends towards a uniform distribution, resulting in 
much lower phase sensitivity. This leads to a significant increase in the incidence of phase 
wrapping, causing the measured vibration amplitude in the phase-sensitive method to tend 
towards zero, thus, underestimating the true value. However, the amplitude PDF continues to 
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be approximated by an unbiased Gaussian distribution with constant variance whenever the 
OCT SNR > -5 dB - a clear advantage of STdOCE. For example, in Fig. 4, at OCT SNR = 5 
dB the mean vibration amplitude measured using the phase-sensitive and STd methods are 
-0.34 urn and -0.42 urn, respectively. This underestimation of the mean vibration amplitude 
using the phase -sensitive method is consistent with previously presented results in flow 
imaging [29]. Using the vibration amplitude measured with the phase-sensitive method at 
OCT SNR = 55 dB as the best estimate of the mirror's true vibration amplitude, at OCT SNR 
= 5 dB, the phase-sensitive method underestimates the vibration amplitude by >20%, whereas, 
STdOCE underestimates it by <3%. 
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Fig. 4. Vibration amplitude (±1 SD) of a mirror measured using STdOCE (blue) and phase- 
sensitive OCE (red) versus OCT SNR. 

It should be noted that for OCT SNR < 5 dB, STdOCE will also introduce significant 
errors in the vibration amplitude measurement as, in this case, the amplitude PDF will 
approach a Rayleigh distribution. The pertinent point is that this error occurs for much lower 
SNR than in the phase-sensitive method. 

For OCE imaging of turbid samples, in which the OCT SNR decreases with depth due to 
attenuation, the results in Fig. 4 indicate that STdOCE should provide accurate measurement 
of vibration amplitude over a larger depth range than the phase-sensitive method, which we 
shall confirm in the next section. This result also indicates that in elastograms generated using 
the phase-sensitive method, the strain calculated in regions of low OCT SNR (<20 dB) will be 
artificially high, due to underestimation of the vibration amplitude. This is discussed in more 
detail in the next section. 

4.2. Elastograms 

In this section, we present 2D elastograms generated using both STd and phase-sensitive 
methods applied to the same data. In each case, the vibration amplitude of the unloaded ring 
actuator was measured using both methods to be 0.35 urn. 

4.2.1. Phantom 2: Soft with hard inclusions 

A structural OCT image of Phantom 2 is shown in Fig. 5(a). The inclusion is located in the 
center of the image and is indicated by the labeled arrow in the figure. Below the inclusion, a 
shadow artifact, corresponding to a region of low OCT SNR, is also labeled. A plot of the 
OCT SNR at the lateral position indicated by the vertical red arrow in Fig. 5(a) is shown in 
Fig. 5(d). In this plot, the large slope of the OCT SNR with depth in the inclusion indicates 
high attenuation, as expected due to the increased concentration of scatterers in this feature. 

Vibration amplitude images generated using the STd and phase-sensitive methods are 
shown in Figs. 5(b) and 5(f), respectively. In both images, the stiff inclusions are visible due 
to a lower rate of change in vibration amplitude with depth, i.e., lower strain, in this feature. 
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Fig. 5. Phantom 2: (a) OCT structural image; (b) Vibration amplitude image; and (c) 
Elastogram for STdOCE; (d) OCT A-scan; and (e) Vibration amplitude plots for STdOCE 
(blue) and phase-sensitive OCE (red) at the lateral position indicated by the red arrow in (a), 
where the dashed lines in (d) and (e) indicate the boundaries between the soft bulk and hard 
inclusion; (f) Vibration amplitude image; and (g) Elastogram for phase-sensitive OCE. 

Vibration amplitude plots generated using both methods, at the lateral position indicated by 
the red arrow in Fig. 5(a), are shown in Fig. 5(e). Both plots match well until a depth of -600 
urn. At depths >600 urn, the decrease in vibration amplitude is higher for the phase-sensitive 
method (red). This is caused by an underestimation of the vibration amplitude in the low OCT 
SNR region below the inclusion (see Subsection 4.1) and results in an artificially high value 
of calculated strain. This result is discussed in greater detail in Section 5. 

The elastograms corresponding to Figs. 5(b) and 5(f) are shown in Figs. 5(c) and 5(g), 
respectively. The artificially high strain in the phase-sensitive elastogram (Fig. 5(g)) is clearly 
visible below the inclusion. As the elastogram is used as a surrogate for elasticity, this leads to 
errors in the interpretation of elastograms. For example, in Fig. 5(g), the region below the 
inclusion is incorrectly presented as having a significantly lower stiffness than the 
surrounding soft material. 

4.2.2. Porcine tissue 

To demonstrate STdOCE on tissue, we performed measurements on freshly excised porcine 
tissue. A flat 10 mm x 10 mm x 2.5 mm section of tissue taken from the chest wall was 
positioned in our OCE system. A representative histology slice taken after OCT/OCE imaging 
is shown in Fig. 6(a). Two tissue types can clearly be identified in this image. A thin soft layer 
approximately 150 urn thick is visible near the surface of the tissue. This layer corresponds to 
adventitia. The underlying stiffer thick layer is intercostal muscle. The corresponding OCT 
structural image is shown in Fig. 6(d). The two layers distinguishable in the histology slice 
can readily be seen in this image: with the thin adventitia layer having a higher OCT SNR 
than the underlying muscle. The vibration amplitude images obtained using STdOCE (Fig. 
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6(b)) and phase-sensitive OCE (Fig. 6(e)) are also shown. The corresponding elastograms are 
shown in Figs. 6(c) and 6(f). The two layers visible in the histology and OCT images can be 
distinguished in both elastograms. Higher strain is present in the softer thin layer in 
comparison to the stiffer underlying tissue, as expected. However, the strain contrast between 
the layers is significantly higher for STdOCE than for phase -sensitive OCE. The strain 
measured in the thin layer is similar using both techniques. However, much higher strain is 
visible in the thicker layer using the phase-sensitive method. This higher strain is caused by 
underestimation of the vibration amplitude measurements using this method due to low OCT 
SNR in the thick layer. The elastograms presented in Fig. 6 demonstrate the improvement in 
quality obtained using STdOCE in comparison to phase-sensitive OCE. 

This advantage can also be seen in Fig. 7, which shows the OCT SNR (Fig. 7(a)) and 
vibration amplitude (Fig. 7(b)) plots corresponding to the lateral position indicated by the 
arrow in Fig. 6(d). Figure 7(a) confirms the high OCT SNR in the thin layer, as well as the 
significant decrease in OCT SNR in the thick layer. Figure 7(b) shows the corresponding 
vibration amplitude plots obtained using STdOCE (blue) and phase-sensitive OCE (red). The 
vibration amplitude measured using both methods matches well until a depth of -480 urn. At 
depths >480 urn, a significantly larger slope (strain) is observed in the phase-sensitive plot. 
This change in slope does not correspond to any structural changes visible in either the 
histology or the OCT images and is caused by the low OCT SNR at these depths. It is this 
artificial increase in strain that reduces the contrast between the two layers in the 
corresponding phase-sensitive elastogram (Fig. 6(f)). In contrast, it can be seen in Fig. 7(b) 
that the STdOCE vibration amplitude plot (blue) maintains a relatively constant slope at 
depths up to 720 urn, indicating that the thick layer is mechanically uniform. For the result 
presented in Fig. 7, we observe a -50% increase in the depth range over which reliable 
vibration amplitude measurements are obtained using STdOCE in comparison to phase- 
sensitive OCE. 




Fig. 6. Porcine tissue: (a) Histology; (b) Vibration amplitude image; and (c) Elastogram for 
STdOCE; (d) OCT structural image. The red arrow indicates the location of the plots in Fig. 7; 
(e) Vibration amplitude image; and (f) Elastogram for phase-sensitive OCE. 
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Fig. 7. (a) OCT A-scan; and (b) Vibration amplitude plots of porcine tissue for STdOCE (blue) 
and phase-sensitive OCE (red), at lateral position indicated by the red arrow in Fig. 6(d). The 
dashed lines in (a) and (b) indicate the boundary between tissue layers. 

5. Discussion 

In this paper, we have presented a new method to extend the performance of dynamic OCE in 
regions of the elastogram corresponding to low OCT SNR. From the experimental data 
presented in Fig. 4, the accuracy of STdOCE exceeds that of phase-sensitive OCE for OCT 
SNRs below ~20 dB. Since the OCT SNR is often below 20 dB, this improvement will lead to 
a marked general improvement in elastograms. 

In Subsection 4.2.1, we experimentally confirmed the systematic underestimation of 
vibration amplitude in regions in a phantom corresponding to low OCT SNR (see Fig. 5). The 
soft material in the inclusion phantom has uniform elasticity [5], thus, assuming uniform 
stress in this material, the strain above and below the inclusion should ideally be the same. 
From Fig. 5(e) (blue curve), the strain measured with STdOCE approximates this ideal case 
well, with strains of 4.1 x 10~ 4 and 3.9 x 10~ 4 measured above and below the inclusion, 
respectively, versus 3.8 x 10~ 4 and 6.1 x 10~ 4 with the phase-sensitive method (red curve). 
Secondly, as the phantom has uniform elasticity at depths beneath the inclusion, i.e., >600 
urn, the true vibration amplitude will decay linearly with depth from this point, reducing to 
zero at the opposite surface of the phantom, i.e., at a depth of 2.5 mm, corresponding to the 
thickness of the loaded sample. If we extrapolate the linear decay in vibration amplitude 
visible in Fig. 5(e) for depths >600 um, we predict, based on STdOCE, zero vibration 
amplitude at a depth of ~2.5 mm, as expected. Performing the same analysis for phase- 
sensitive OCE, we incorrectly predict zero vibration amplitude at a depth of ~2.1 mm, 0.4 mm 
from the sample boundary. 

As comparisons made between phase-sensitive and STd methods are based on the same 
data, the data acquisition time was the same for both methods. STdOCE relies on the analysis 
of multiple frequency tones, thus, dense sampling is required to adequately sample each tone. 
In comparison, phase-sensitive OCE analyzes the phase difference spectrum, consisting of a 
single frequency component. The vibration amplitude can be estimated using the phase- 
sensitive method once the A-scan acquisition frequency exceeds twice the excitation 
frequency [17], at the expense of lower accuracy for less densely sampled data [21]. As the 
acquisition speed of OCT systems increases and more sophisticated signal processing 
hardware, e.g., graphics processing units (GPUs) become available, the dense sampling 
required for STdOCE will become more practical and efficient. 

Many related OCT techniques and applications rely on vibration amplitude measurement: 
for imaging vibrations of the vocal chords [41,42] and the ear [21,43,44]; in functional 
imaging techniques such as magneto-motive OCT [19]; and for photothermal OCT [20]. The 
method presented in this paper could readily be applied in such techniques to improve the 
accuracy of vibration amplitude measurement. 
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6. Conclusions 



We have presented a new method to measure vibration amplitude in dynamic optical 
coherence elastography based on analysis of the Doppler spectrum. We have demonstrated 
that our method is more robust for low OCT SNR than the currently used phase-sensitive 
method. Vibration amplitude measurements performed on a vibrating mirror when the OCT 
SNR was held constant at 5 dB demonstrate <3% error in comparison to >20% using the 
currently employed phase-sensitive method. Elastograms of a tissue-mimicking phantom and 
excised porcine tissue demonstrate improvements in quality with this new method, including a 
50% increase in the depth range of reliable vibration amplitude. Such improvements increase 
the clinical potential of this emerging modality. 
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